INTERFACE FOLIATION NEAR MINIMAL SUBMANIFOLDS IN 
RIEMANNIAN MANIFOLDS WITH POSITIVE RICCI 

CURVATURE 



MANUEL DEL PINO, MICHAL KOWALCZYK, JUNCHENG WEI, AND JUN YANG 

Abstract. Let {M,g) be an A'^-dimensional smooth compact Riemannian 
manifold. We consider the singularly perturbed AUen-Cahn equation 

e^AgU + (1 - u'^)u = in A4, 

where e is a small parameter. Let K C M be an (A'^ — l)-dimcnsional smooth 
minimal submanifold that separates Ai into two disjoint components. Assume 
that K is non-degenerate in the sense that it docs not support non-trivial 
Jacobi fields, and that |^/cP + RiCg(!^Ki i'k) is positive along IC. Then for 
each integer m > 2, we establish the existence of a sequence e = ej 0, and 
solutions ttg with m-transition layers near K, with mutual distance 0(e| Ine]). 
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1. Introduction 



In the gradient theory of phase transitions by Allen-Cahn [2], two phases of 
a material, +1 and —1 coexist in a region il C separated by an {N — 1)- 
dimensional interface. The phase is ideahzed as a smooth e-regularization of the 
discrete function, which is selected as a critical point of the energy 



where e > is a small parameter. While any function with values ±1 minimizes 
exactly the second term, the presence of the gradient term conveys a balance in 
which the interface is selected asymptotically as stationary for perimeter. The 
energy may be regarded as an e -relaxation of the surface area: indeed, in [25] 
it is established that a sequence of local minimizers Wg, with uniformly bounded 
energy, must converge in Ll^^-sense to a function of the form xe — XE" so that 
dE locally minimizes perimeter, thus being a (generalized) minimal surface. This 
is the starting point of the T-convergence theory, in which the constraint of to 
a suitable class of separating-phase functions, converges to the perimeter function 
of the interface. Indeed, analogous assertions hold true for general families of 
critical points, and for stronger notions of interface convergence, see [6, 29, 33]. The 
principle above applies to modeling phase transition phenomena in many contexts: 
material science, superconductivity, population dynamics and biological pattern 
formation, see for instance [31] and references therein. 

It is natural to consider situations in which phase transitions take place in a 
manifold rather than in a subset of Euclidean space. In this paper we consider 
a compact A^-dimensional Riemannian manifold {M , g) , and want to investigate 
critical points in H^{Ai) of the functional 



with sharp transitions between —1 and 1 taking place near a {N — l)-dimensional 
minimal submanifolds of A4. Critical points of correspond precisely to classical 
solutions of the Allen-Cahn equation in Ai, 



where Ag is the Laplace-Beltrami operator on A4. 

Wc let in what follows /C be a minimal (A^ — l)-dimcnsional embedded submani- 
fold of Al, which divides Ai into two open components Ai±. (The latter condition 
is not needed in some cases.) The Jacobi operator J of /C, corresponds to the second 
variation of A^- volume along normal perturbations of /C inside A\ : given any smooth 
small function v on /C, let us consider the manifold IC{v), the normal graph on /C 
of the function v, namely the image of JC by the map p £ K, i-^ expp(v{p)i']c{p)) . 
If H{v) denotes the mean curvature of lC{v), defined as the arithmetic mean of 
the principal curvatures, then the linear operator is the differential of the map 
V I—!- nH{v) at w = 0. More explicitly, it can be shown that 






(1.1) 



(1.2) 
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where A^; is the Laplace-Beltrami operator on /C, \A]c\'^ denotes the norm of the 
second fundamental form of K., RiCg is the Ricci tensor of M and vjc is a unit 
normal to /C. We will briefly review these concepts in Section 2. 

The minimal submanifold JC is said to be nondegenerate if the are no nontrivial 
smooth solutions to the homogeneous problem 

J^jj = in/C. (1.3) 

This condition implies that K. is isolated as a minimal submanifold of A4. 

In [28], Pacard and Ritore assume that /C is non-degenerate and, and proved 
that there exists a solution to equation (1.1) with values close to ±1 inside 
Ai±, whose (sharp) 0-level set is a smooth manifold which lies e-close to /C. More 

precisely, let w{z) := tanh (^-^^ be the unique solution of the problem 

w" + w-w^ = in M, w{0) = 0, w{±oo) = ±1, (1.4) 
and denote by c* its total energy, namely 

Then the solution in [28] resembles near K, the function w{t/£), where f is a 
choice of signed geodesic distance to F. In particular 

Je{Ue) C* |/C| . 

In this paper we describe a new phenomenon induced by the presence of positive 
curvature in the ambient manifold M: in addition to non-degeneracy of /C, let us 
assume that 

K - \Ak.\'' +mCg{v^,VK.)>^ on/C. (1.5) 

Then, besides the solution by Pacard and Ritore, there are solutions with multiple 
interfaces collapsing onto K,. In fact, given any integer m > 2. we find a solution 
Ue such that — 1 approaches in ^A± as e — > 0, with zero level set constituted 
by m smooth components 0(e| log£|) distant one to each other and to /C, and such 
that 

Je{Ue) mc* \K.\ . 

Condition (1.5) is satisfied automatically if the manifold M has non-negative Ricci 
curvature. If iV = 2, _ftr corresponds simply to the Gauss curvature of Ai measured 
along the geodesic JC. 

The nature of these solutions is drastically different from the single-interface 
solution by Pacard and Ritore [28]. They are actually defined only if e satisfies a 
nonresonance condition in s. In fact, in the construction e must remain suitably 
away from certain values where a shift in Morse index occurs. We expect that 
the solutions we find have a Morse index 0(|loge|°) for some a > as critical 
points of Je, while the single interface solution is likely to have its Morse index 
uniformly bounded by the index of K. (namely the number of negative eigenvalues 
of the operator J). 
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Theorem 1. Assume that K, is nondegenerate and embedded, and that condition 
[1.5) is satisfied. Then, for each m> 2, there exists a sequence of values e = £j — > 
such that problem {1.1) has a solution such that w^ — 1 — >■ uniformly on compact 
subsets of M±, while nearK,, we have 



1=1 



I) 



1 



where (jj, z) are the Fermi coordinates defined near K. through the exponential map 
(see Section 2.1), and the functions ft satisfy 



o(i) 



(1.6) 



Moreover, when N ~ 2, there exist positive numbers i>i, 
c > and all sufficiently small e > satisfying 

I Vi c 

> — , for all 



log- 



J' 



1, 



J" 



, , m — 1, 



, Vm-i such that given 
= 1,2,.... (1.7) 



a solution Ue with the above properties exists. 



We observe that the same result holds if m is even and M\K. consists of just one 
component. Thus the condition that /C divides M. into two connected components 
is not essential in general. 

As we will see in the course of the proof, the equilibrium location of the interfaces 
is asymptotically governed by a small perturbation of the Jacobi-Toda system 

e^(^A,cfj + {\AK\^ + RiCg{i^K,'yK))fj) - ao[e'^^'"^^-'^ - e-^f'+'~f^^] = 0, (1.8) 

on K., j ~ l,...,m, with the conventions /o = — oo, /,n+i = +oo. Heuristically, 
the interface foliation near /C is possible due to a balance between the interfacial 
energy, which decreases as the interfaces approach each other, and the fact that the 
length or area of each individual interface increases as the interface is closer to K, 
since M is positively curved near K.. 

What is unexpected, is the need of a nonresonance condition in order to solve the 
Jacobi-Toda system. A question which is of independent interest is the solvability 
of the Jacobi-Toda system without the condition (1.5). Similar resonance has been 
observed the problem of building foliations of a neighborhood of a geodesic by CMC 
tubes considered in [17, 22]. This has also been the case for (simple) concentration 
phenomena for various elliptic problems, see [8, 16, 19, 20]. 

Our result deals with situations in which the minimal submanifold is locally but 
not globally area minimizing. In fact, since condition (1.5) holds, the Jacobi op- 
erator has at least one negative eigenvalue, and near /C, cannot have parabolic 
points. In the case of a bounded domain Q. of under Neumann boundary condi- 
tions, a multiple-layer solution near a non-minimizing straight segment orthogonal 
to the boundary was built in [9]. In ODE cases for the AUen-Cahn equation, 
clustering interfaces had been previously observed in [7, 26, 27]. No resonance phe- 
nomenon is present in those situations, constituting a major qualitative difference 
with the current setting. 
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The method consists of hneariziiig the equation around the approximation 

"o(x,.) = f:-(^^f^) + ^((-ir-^-i). 

and then consider a projected form of the equation which can be solved boundedly 
after finding a satisfactory hnear theory, and then applying the contraction mapping 
principe. In that process the functions fj are left as arbitrary functions under some 
growth constraints. At the last step one gets an equation which can be described 
as a small perturbation of the Jacobi-Toda system 

We do not expect that interface foliation occurs if the limiting interface is a 
minimizer of the perimeter since in such a case both perimeter of the interfaces and 
their interactions decrease the energy, so no balance for their equilibrium locations 
is possible. On the other hand, negative Gauss curvature seems also prevent inter- 
face foliation. This is suggested by a version of De Giorgi-Gibbons conjecture for 
problem (1.1) with Ai the hyperbolic space, established in [3]. 

2. Geometric background and the ansatz 

In the first preliminary part of this section, we list some necessary notions from 
differential geometry: Fermi coordinates near a submanifold of A4, minimal sub- 
manifold, as well as Laplace-Beltrami and Jacobi operators. We then express the 
problem in a suitable form, define an approximate solution and estimate its error. 

2.1. Local coordinates. Let be an iV > 2-dimensional smooth compact Rie- 
mannian manifold without boundary with given metric g. We assume that K. is 
an iV — 1 dimensional submanifold of A4. For each given point p E K,, TpAi splits 
naturally as 

TpM = TpJC ® NpIC, 

where TpK. is the tangent space to K, and NpK, is its normal complement, which 
spanned respectively by orthonormal bases {Ei : z = 1, • • • , N— 1} and {En}. More 
generally, we have for the tangent and normal bundles over IC the decomposition 

TM = T/C ffi NIC. 

Let us denote by y the connection induced by the metric g and by sj^ the corre- 
sponding normal connection on the normal bundle. 

Notation: Up to section 2.4, we shall always use the following convention for 
the indices 

i,j,kj---e {l,2,--- ,N-1}, a,b,c,---e {1,2,- ■■ ,N}. 

Given p € IC, we use some geodesic coordinates y centered at p. More precisely, 
in a neighborhood oi p in K., we consider normal geodesic coordinates 

y = Yp{Y) ■.= exp^{Y,E,), y = (yi, • • • , y^^i) G V, (2.1) 

where exp'^ is the exponential map on IC and summation over repeated indices is 
understood. V is a neighborhood of the origin in R^^^. 

This yields the coordinate vector fields Xi ~ f*{di), i ~ 1, ■ ■ ■ , N — 1 where 
/(y) = Yp{f). For any E € Tp/C, the curve 

s Je{s) = cxpp{sE) 
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is a geodesic in /C, so that 

Vx,X.j \p e NpIC for any i, j = 1, • • • , iV - 1. 

We recall that the Christoffel symbols FjJ, i, j — 1, ■ ■ ■ ,N~1 arc given by 

Vx,X,\p ^ T^jEn, i. c. r,^ ^ ~g{Vx,X,,EN). 

We also assume that at p the normal vector E^ is transported parallelly (with 
respect to V^) through geodesies 7b (s) from p. This yields a frame field for 
NIC in a neighborhood of p which satisfies 

Vx,^jv|p e Tp/C, i.e. ~g{yE,EN,EM)\p = 0, i = l,---,N-l. 

We define the numbers F^j^' *7.7 = 1; ' ' ' j ^ — Ij by 

Af-l 

Vx.Xjvlp = Y^^'^nE,, i.e. = g(Vx,Xjv,i?,). 

In a neighborhood of p in A^, we choose the Fermi coordinates (y, 5) on Ai 
defined by 

$°(y,z) = expy^(y)(zi;Ar) with (y, 5) = (yi,-" ,yw-i,^) e V x {-5o,5a), (2.2) 

where expy^^y-j is the exponential map at Yp{y) in Ai. We also have corresponding 
coordinate vector fields 

X, = $^(%), Xx^^lidi). 
By construction, X^rlp = E^- 

2.2. Taylor expansion of the metric. In this section we will follow the notation 
and calculations of [17]. By our choice of coordinates and the Gauss Lemma, on K. 
the metric g splits in the following way, 

Af-l 

9{p) = X! dij'^yi ® + gNNdz (g) dz, p e /C. (2.3) 

As usual, the Fermi coordinates above are chosen so that the metric coeflicients 
satisfy 

9ab = g{Xa, Xb) = 6ab at p. 

Furthermore, g{Xi, X^) = in some neighborhood of p in /C. Then 

Xig{Xj, Xn) = g{\/x,Xj, Xn) + g{Xj, \/x.Xn) on /C, 
yield the identity 

riI + r^A. = o at p. (2.4) 

We denote by F^ : NJC — > M, = 1, • • • — 1, the 1-forms defined on the 
normal bundle of JC as 

TI{En) = g{VE,E,,EN). (2.5) 

The second fundamental form Ajc : TIC x T/C — > A^/C of the submanifold K- and its 
corresponding norm are then given by 

AT-l 2 

A^{E,,E,) = T\{Ej,)Ej,, = Yl ■ (2.6) 
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For X, Y,Z,We TM , the curvature operator and curvature tensor are respectively 
defined by the relations 

R{X, Y, Z)^^xVy Z-^yVxZ - ^[x,y]Z, (2.7) 
R{X,Y,Z,W)^g{R{Z,W)Y,X). (2.8) 

The Ricci tensor oi {Ai,g) is defined by 

mCg{X,Y)^r''R{X,Xa,Y,Xt,). (2.9) 

We now compute higher order terms in the Taylor expansions of the metric 
coefficients. The metric coefficients at g = $'^(0, z) are given in terms of geometric 
data at p = 4>°(0, 0) and \z\ = distg(p, g), which is expressed by the next lemmas, 
see Proposition 2.1 in [17] and the references therein. 

Lemma 2.1. At the point q = $'^(0, z), the following expansions hold 

Vx„ Xn = Oi\S\)Xa, (2.10) 





Vx^Xj=Tl{EN)XN + 0{\z\)Xa, ^,J = 1,••• ,A^-1, (2.11) 



Af-l 

Vx,XN^Vx.X, = Y,^iiEN)X,+0{\S\)Xa, i = l,--- ,iV-l. (2.12) 

i=i 

Lemma 2.2. In the above coordinates {y, z), for any i,j = l,2,-- - ,A^— 1, we have 



5.y(0,z) = S,,-2rl{EN)S~RiXN,X,,XN,X,)\2 



z\2 



J2rHEN)n{Ex)\S\' + 0{\z\'), (2.13) 



k=l 

5.jv(0,z) = Od^n, (2.14) 

5wjv(0,z) = l + 0(|z|3). (2.15) 

2.3. The Laplace-Beltrami and Jacobi operators. If {Ai,g) is an A^-dimcnsional 
Riemannian manifold, the Laplace-Beltrami operator on Ai is defined in local co- 
ordinates by the formula 

1 



Am = ^^da[Vd^r''db), (2.16) 

where g"'^ denotes the inverse of the matrix (gab)- Let /C C be an (N — 1)- 
dimensional closed smooth embedded submanifold associated with the metric go 
induced from {A4,g). Let A^; be the Laplace-Beltrami operator defined on /C. 

Let us consider the space C°°{N1C), which identifies with that of all smooth 
normal vector fields on /C. Since /C is a submanifold of codimension 1, then given 
a choice of orientation and unit normal vector field along /C, denoted by vjc G NIC, 
we can write e C°°{N]C) as = 01/^, where G C°°{JC). 

For G C°°{NJC), we consider the one-parameter family of submanifolds t ^■ 
/Ct^* given by 

ICt,^ = {cxpy{t^iy)):yeJCy (2.17) 
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The first variation formula of the volume functional is defined as 

^ Vol(/Ct.*) - / < *, h >w dVK, (2.18) 
dt t=o Jk. 

where h is the mean curvature vector of /C in , < • , • > jv denotes the restriction 
of g to NIC, and AVk. the volume element of K,. 

The submanifold K, is said to be minimal if it is stationary point for the volume 
functional, namely if 
d 



dt 



Vol(/Ctvi,) = for any * e C°°(iV/C), (2.19) 

t=o 

or cquivalently by (2.18), if the mean curvature h is identically zero on JC. It is a 
standard fact that if TI{Em) is as in (2.5), then 

K is minimal <^=^ ^ ^\{En) = 0. (2.20) 

i=l 

The Jacobi operator J appears in the expression of the second variation of the 
volume functional for a minimal submanifold K, 

d2 



dt2 



t=0 



Vol{ICt.^) = - J < J^, ^ >N dVic for any * e C°°(A^/C), (2.21) 
and is given by 

= ^A^(l) - RiCgivK, i^k)(I) - \AK\^(j), (2.22) 

where ^' ~ (p^Ki as has been explained above. 

The submanifold K, is said to be non-degenerate if the Jacobi operator is 
invertible, or equivalently if the equation J'(j) = has only the trivial solution in 
C°°{JC) 

2.4. Laplace-Beltrami Operator in Stretched Fermi Coordinates. To con- 
struct the approximation to a solution of (1.1), which concentrates near /C, after 
rescaling, in AA/e, we introduce stretched Fermi coordinates in the neighborhood 
of the point e~^p € £~^1C by 

$,(y,z) = i$0(ey,ez), (y, z) = (y^, . . . , y^.^, z) e g-iy x ( - ^, ^) . (2.23) 

Obviously, in Ai^ = e~^A4 the new coefficients gat's of the Riemannian metric, 
after rescaling, can be written as 

9ab{Y, z) = 5afc(ey, a, & = 1, 2, • • • ,N. 

Lemma 2.3. In the above coordinates (y, z), for any i, j — 1,2, ■ ■ ■ , N —1, we have 

g,,{j,z) = 51 ^2eT{(EN)z - e^R(Xj,XN,XN,Xi)\z\'' 

N-l 

+ e^Y. rf (i?iv)ri.(i?^)|z|2 + 0{\ez\^), (2.24) 



k=l 

|2^ 



g,N{y,z) ^ 0{\ez\'), (2.25) 
gNN{y,z) = l + 0{\ezf). (2.26) 
Here R{-) and F^ are computed at the point p ^ JC parameterized by (0,0). 
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Now we will focus on the expansion of the Laplace-Beltrami operator defined by 

Vdetg V / ^227) 

= g''''dadb + {da9'''')dt + ia4log(det5))5°'afc. 

Using the assumption that the submanifold K. is minimal as in formula (2.20), direct 
computation gives that 

detg = 1 - e^K{ey)z^ + 0{e^\z\^), 
where we have, using (2.6) and (2.9), denoted 

K = mCg{v^,VK,) + \AK.\^. (2.28) 

This gives 

log (detg) = -e^K{ey)z^ + 0{eM'')- 
Hence, we have the expansion 

Am, = 9.. + A,c, + e'^zKiey) + B (2.29) 
where the operator B has the form 

B ^ ez 5y + e^z^ af^ d,, + e^z^ a%j^ d,, + e^z bid., + e^z^ b% d, . (2.30) 

and all the coefficients are smooth functions defined on a neighborhood of K, in , 
evaluated at {ey,ez). 

2.5. The local approximate solution. If we set u{x) := u{ex), then problem 
(1.1) is thus equivalent to 

Am, u + F{u) =0 in M^, (2.31) 

where F{u) = u — . In the sequel, we denote by and /C^ the e~^-dilations 
of M and JC. 

To define the approximate solution we observe the heteroclinic solution to (1.4) 
has the asymptotic properties 

«;(z) - 1= -2e-^l^l + 0(e-2V2|^l), z > 1, 

w;(z) + 1= 2e-^l^l + 0(e-2^l^l), z < -1, (2.32) 

w'{z) = 2V2e-^l^l + 0(e-2^l^l), \z\ > 1, 

For a fixed integer m > 2, we assume that the location of the m phase transition 
layers are characterized in the coordinate {y, z) defined in (2.23) by the functions 

^ = /j (£y)7 1 < j < with 

/i(£y) < f2{£y) < ■ ■ ■ < frniey), 

separated one to each other by large distances 0{ \ loge|), and define in coordinates 
(y, z) the approximation 

uo{y,z) ■.= Y^w,{z^ f,(ey)) + , w,{t) := {-iy-^w{t), (2.33) 

with this definition we have that uo(2/, z) Wj{z — fj{ey)) for values of z close to 
fjiey). 
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The functions fj : /C — > M will be left as parameters, on which we will assume a 
set of constraints that we describe next. 

Let us fix numbers p > TV, M > 0, and consider functions hj E W^'P{IC), 
j = 1, . . . , m, such that 

\\hj\\w^.p{K) ■= II-D^^jIIlpCC) + \\Dichj\\LP{K) + < M . (2.34) 

For a small e > 0, we consider the unique number p = with 

e-^p = e^p. (2.35) 
We observe that p^ is a large number that can be expanded in e as 



£ \/2 V ^eJ wog _ 
Then we assume that the m functions /j : /C — > M arc given by the relations 

/ TO + 1 \ 

fkiy) = ik —\p^ + hk{y), k = l,...,m, (2.36) 

so that 

fk+i{y) - fkiij) = Pe + hk+iiy) - hkiij) , fc = l,2,...,m-l. (2.37) 

We will use in addition the conventions ho = — oo, hm+i = +oo. 

Our first goal is to compute the error of approximation in a Sq j s neighborhood 
of /Ce, namely the quantity: 

S(un) = Am^uo + i^(wo). (2.38) 

For each fixed 1 < ^ < to, this error reproduces a similar pattern on each set 
of the form 

M = ^ (y, z) e /C, X (-^, / k - !i[ey)\ < ip, + A/ I . (2.39) 

For (y, z) G Ag^ we write t = z — fi{£y) and estimate in this range the quantity 
S{uo)(y,t + f({ey)). We have the validity of the following expression. 

Lemma 2.4. For £ E {1, . . . , m} and {y, z) € A( we have 
{-lY-'S{uo){y,t + ff) 

= 6{l -w\t))e^pe re-^('*''-''*-i)e^* - e'^^'^'+'^'^'^e-^'] (2.40) 
^e' (^A,che + it + h)K)w\t) + e^\^K h\''w"{t) + {-lY-^ee{ey,t) . 
where for some r, cr > we have 
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Proof. From (2.29), using that w'- + F{wj) — 0, we derive that, for {y, z) e Ai 

m m 

S{uo) =F{m{y,z)) - ^F(i«,(z-/,(ey)) + ^ | /i,(£y)|'<(^ - /.N)) 
j=i i=i 

-e'Y. i^^'hi^y) + ^K{ey) )w'^{z - f,{ey)) 



i=i 

+ [a\k{^y,£z)d^khj{£y) + hl{ez,ey) dihj{ey)] Wj{z - fjiey)) 

+ e^[z^ a%j^{ey,ez) + z'^ afj^{ey,ez) d^hj{ey)] (z ~ fj{ey)) 

+ z al^. {ey, ez) dih, (ey) dkhj [ey) w'- {z - fj (ey)) . (2.41) 
Lsider first the cas 
We begin with the term 



3 

Let us consider first the case 2 < ^ < m — 1 . 



F{uo{y,t + fe)) - Y.F{w,{t + fe - f,)) , \t\ < ^. 

Since 

w{s) = ±(1 - 2e-^l^l) + 0(e-2^l"l) as s ±oo, 
we find that for j < £, 

w{t + h- fj) - 1 = ~2e-^^f'-f^^e-^' + O (e-2y2|t+/,-/,| j ^ ^2.42) 
while for j > i, 

w{t + fe - fj) + 1 = 2e-^^f'-f'^e^' + O (e-2V2|t+/.-/,| ^ _ (2.43) 
Now, since 

F{w)^w{l^w^), |t|< 1 + 0(1), 

\fi-fi\ = \t^]\Pe + 0{\), e-^^P' ^e^p„ 
wc find that if \j ~ (\ > 2 and < a < \/2, then for sonic r > 0, 

\F{wj{t + fi- f,))\ < C7e-^l*+^*"^^l < e^+^e-'^l*!. (2.44) 
On the other hand, for certain numbers Si, S2 G (0, 1) we have 

F{w{t + f( - /,_i)) = F'{1) ai, +iF"(l + siai) aj, (2.45) 

+ - /£+i)) = i^'(l) 02 + ^^^"(1 - 8202) 4. (2.46) 

where 

ai := !«(< + - /f-i) - 1, a2 := w(i + fi - fi+i) + 1- 

Now, we find 



(—1)^ ^un = w(t) — fli — a2 — 03, 03 = O f max e 

Vb-«l>2 



V^|t+/f-/j| 
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Thus for some S3 G (0, 1), 



(-l)^-iFK) ^Fiw) - F'{w)iai+a2) + - 53(01 + a2))(ai + 02)' 



+ O ( max e-^l*+^'"^^l 
vb-^l>2 



Combining relations (2.44)-(2.47) and using that 

F'{1) - F'iw) = 3(1 - w'), |ai| + \a2\ - 0(e-^*) ^ 0(e-^l*l), 
we obtain 



(2.47) 



= 3(1 (01+02) 



F"{1 - 5202) - F"{w ~ .S3(ai + a2))] (a? + 0^ + (e'+^e"'^!*! ) 



3(l-u;2)(aj+a2) + O (e^+^g— 1*1 ) 



Hence, recalling relations (2.37), (2.35), the definitions of ai, 02 and the asymptotic 
expansions (2.42), (2.43) for j ~ £ — 1 and j = ^ + 1, we find 



F(uo)-^F(«;,(t-/,)) 



6(~l)^-i(l~i.2(t))e2p, 



where = O (e^+^e"'"!*! ) . 

Substituting (2.48) in expression (2.41) wc then find 

{-lY-^Siuo){y,t + fe) 
= 6{1 ~w^t))e^p, 

- (^A,che + (t + fe)K)w\t) + e^l^ic he\^w'{t) 
Here we have denoted 



g-\/2(/i,-?i,_i)g\/2t _ g-\/2(/i«+i-/if)g-y2t 



(2.48) 



(2.49) 
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e,(ey, t) = e,{ey, t)- e^J2 {^i^''^ + + /^)^ ) "'K^ + " 



+ £3(-l)^-iz fli, d,h, dkhj w';{t + fe- fj) 

+ e^it + /,) [a]^ d^khj + bl d,h^ ] w]{t + ,h - fj) 

+ e^[z^ a^jv + oIn ^^hj\ w]it + /, - fj), (2.50) 

where the coefhcients are understood to be evahiated at ey or (ey, £{t + f({ey))). 

While this expression has been obtained assuming 2<£<m — 1, we see that it 
also holds ior £ = m, £ = 1. The cases £ = 1 and £ = m are dealt similarly. The 
only difference is that the term (,~V2{hi-hi-i)^V2t _ ^-V2{he+i-hi)^-V2t gg|-g 

respectively replaced by 

-e-^''^e-^\£=l and e'^'^-e^* , € = m. (2.51) 



2.6. Size of the error. Examining expression (2.50) we see that the error in the 
considered region is made up by terms that can be bounded by a power of e times a 
factor with exponential decay in t. We introduce the following norm for a function 
g{y, t) defined on Kg x R. Let cr > 0, I < p < +oo. We set 

liffllp-- sup e'^l*l||.9||^. y (2.52) 

We want to consider the error associated to points in the set Ai as a function 
defined in the entire space /C^ x M. To do so, we consider a smooth cut-off function 
C(s) with C(s) = 1 for s < 1 and ({s) = for s > 2 and define 

Ce(t)-C(|t|-y-2M). 
We extend the error as follows. Let us set 

Seiuo) := 6(l-u;2(i))eVe [ e-v^(''"-''"-i)eV^* _ g-v^(/.,+i-?.,)g-v^t 1 

+ {-lY-^Ce{t)&e{ey,t), (2.53) 

where the cut-off expressions are understood to be zero outside the support of . 
We see that 

i-lY-^Siuo){ey,t + fe)=Seiuo){y,t) for all (y^t) £ A^. 
The following lemma on the accuracy of the error is readily checked. 
Lemma 2.5. For a given < a < \/2 and any p > 1 we have the estimates 

\\Si{uo)\\.p,a < C£^-\ WCe^iWp.. < Ce^-^ (2.54) 
where r is any number with r > and t > . 
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Proof. The proof amounts to a straightforward verification of the bound term by 
term. Let us consider for instance 

Then bi \t\ < ^ we get 

\Ei\ < Ce^llogel < Ce-^^' < g--^^ e-^^-^"")"- < e"'^l*l e^^^, 

where r > ^^cr- This impHcs llE'ilIp,^ < Ce^^"^ for any 1 < p < +oo. Now, let us 
consider the term 

E2{y,t)^e^ Aichf {ey)w'{t). 
Then for any cr < \/2 wc have 

e'"'''||i^2||Lp(i3((t,y)4) < C!e^\\AK:hi{ey)\\LP{B{yA) 

< Ce^-^\\A^hi\\LPiK) < Ce^-^. 
The rest of the terms are dealt similarly, being in fact roughly at least e times 
smaller than those above. □ 

Very important for subsequent developments is the Lipschitz character of the 
error in the parameter function h = (/ii, . . . , /iat). Let us write Sjih) to emphasize 
the dependence on this function. We have 

Lemma 2.6. Let us assume that the vector-valued functions hi, h2 satisfy the 
constraints (2.34). We have the validity of the following Lipschitz conditions. 

||5,(hi)-5,(h2)||p,, <Ce2--||hi-li2 



-2|1 



for T > and r > 

Proof. Again the proof consists in establishing the bound for each of its individual 
terms, more precisely, we need to bound now for instance dd-h^Sjija.). Since the 
dependence on this object, and as well on second derivatives comes in linear or 
quadratic way, always multiplied by exponentially decaying factors and small pow- 
ers of e, the desired result directly follows. The dependence on the values of the 
functions hj appears in a more nonlinear fashion, however smooth and exponen- 
tially decaying. We omit the details. The complete arguments are rather similar 
to those in the proof of Corollary 5.1 of [10]. □ 

2.7. The global approximation. The approximation uq is so far defined only in 
a neighborhood of K.^ in Me, where the local Fermi coordinates make sense. Let 
us assume that m is an odd number. In this case we require that IC^ separates Me 
into two components that we denote Mj and Mf. 

Let us use the convention that the normal to Ke points in the direction of 
Let us consider the function H defined in Me\ ICe as 

Then our approximation ug (x) approaches M{x) at an exponential rate 0(e~^l*l ) as 
\t\ increases. The global approximation we will use consists simply of interpolating 
Mo with H sufficiently well-inside Me \ ICe through a cut-off in \z\. Let Afg be the 
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set of points in Ms that have Fermi coordinates {y,z) well-defined and \z\ < -. 
with some positive constant S < 6q/10. 

Let ri{s) be a smooth cut-off function with 77(5) = 1 for s < 1 and — for s > 2 
and define 

-<^' lilt 

Then we let our global approximation w(a;) be simply defined as 

V := rjsuo + (l - vs)M, (2.57) 

where EI is given by (2.55) and w is just understood to be EI(a:) outside Afs- 

Since H is an exact solution in \.Me, the global error of approximation is 
simply computed as 

^(w) = Aw + F(w) = rjsS{uo) + E, (2.58) 

where 

E = 2V%Vwo + Avsiua - H) + F{r]sUo + {I - %)H) - %i^K) ■ 

Observe that E has exponential size 0(e~ = ) inside its support, and hence the 
contribution of this error to the entire error is essentially negligible. 

If m is even, we simply define 

^ -.^ VsuQ + {l-vs)i-l). (2.59) 

In this case there is no need that K. separates A4 into two components. 

3. The gluing procedure 

Once the global approximation w(x) in (2.57) or (2.59) has been built, we then 
want to find a solution to the full problem of the form 

u{x) = w(x) -I- ip{x), X G A^e 

where (p{x) is a small function. Thus ip must satisfy 

Am,V + F'{^)v> = -Si^)- N{ip) in Me (3.1) 

where 

N{ip) = F{v + (p)- F(w) - F'{v)ip. 
We shall look for a solution of the form 

rrt 
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where the functions (pj are defined in the entire space /Ce x R. Then the equation 
is equivalent to 



m 

Y,Cj2[^mJj + F'{^)4>, + Oi(F'(w) + 2)V + + + ^(w) 

m rn 
m ??i m 

+ (1 - E Cji) ^(^' + E C«2<^.) + (l - E O2) ^(") = in Me- 

3=1 i=i i=i 

This system wiU be satisfied if the (m + l)-tuplc . . . , ip) solves the system 

AmJj + + Oi(-F'(w) + 2)^ + OiiV(V + <^j) + 5(w) = 0, (3.2) 

for 1^1 < C| loge|, j — I, . . . ,m, and 

- A^^V + 2^^ = Q{i^,x) in TWe, (3.3) 

where we have denoted 

Q(^, x) := (1 - E Oi) 1 [2 + ^' («)] ^ + ^(^ + E ^'2^0 (3-4) 
i=i I 4=1 J 

m m 
+ (1 -E^^-2)5(w) + J2 [2{^MX32,yMj3)+PjAM,Q2 ■ 

The gluing procedure consists in solving equation (3.3) for ?/' in terms of a given 
(j) = (01, . . . , pm) chosen arbitrary but sufficiently small, and then substituting the 
result in equation (3.2). Let us assume the following constraints on the 0j's: 

4>3{y,z) -■ (l>j{y,z~ fj{ey)), WPjh.p.a < I for all j = l,...m,. (3.5) 

Lemma 3.1. Given functions <j)j andh satisfying respectively constraints (3.5) and 
(2.34-), there exists a unique solution tp = ^'(0, h) to equation (3.3) with 

for a small t > 0. In addition the operator ^ satisfies the Lipschitz condition 

||*(0i,hi)-*((/)2,h2)||^ < Ce''--[\\cP^ - cf'h^^^, + - h2h,p]. (3.6) 
Proof. Let us consider first the linear equation 

- /^MA + '2.i^ = E{x) mMe. (3.7) 

We claim that if wc set 

\\E\\p.O = sup |li?|iLP(B(2:,l))i 

then problem (3.7) has, for all small £ > 0, a unique bounded solution ip = 2l(i?), 
which in addition satisfies 

\\DmM + UWoo < C\\E\\p,o, 
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provided that p > m. To prove this claim, it suffices to estabhsh the a priori 
estimate in L°°-norm. If tliat was not true, there would be sequences e = Sm ipm 
En-, with |li?„||p,o 0, llV'nIlcc- = 1 such that 

-Avu^Vn + ZV),! = En in Me- 

Using local normal coordinates around a point p„ £ where |'0n(j'n)l ~ 1' 
same procedure as in the proof of the a priori estimate in Proposition 4.1 leads us 
to local convergence of ipn to a nontrivial bounded solution of 

-AifivV' + = inM^, 

and a contradiction is reached. 

To solve equation (3.3) we write it in fixed point form as 

^-2t(Q(V;,-)). (3.8) 

In the region where the functions (1 — C,ii)C,2j, '^MX'^ji '-^M^C'^j are supported 
we have, thanks to (3.5), 

|0(a^)| + |V^,0(x)| < Ce"^l«^-«-il||0,||2,p,. < e^-^U,h,p,., 
for a small r > 0. We also notice that 

m 

\Q{0,x)\ < Ce^-^ {\Dlh{ey)\ + |i?^h(ey)| + \h{ey)\ + l) + e^-^ ^ |10,|l2,p,.. 

i=i 

We observe then that 

We check next the Lipschitz character of this operator, not just in ip , but also in 
the rest of its arguments. Let us write Q = Q{tp,'h, (p) and assume 

m2,p^a<l, UUp.. < lie^-\ Mw^-^(K)<M. (3.9) 

We consider ('0',0',h'), I = 1,2, satisfying (3.9), and denote Q' = 2(1/)', 0', h'). We 
will show that 

\\Q{ij\<t>\^')-Q{ij\<t>^Y.^)\Wp 

< Ce^-^[U'l - i'2\\^ + - (l)^\\2,p.a + Will - hi\\2,p]. (3.10) 

Let us observe that for (■;/',(/), h) satisfying (3.9), 

Q{x) = Q{i^,(l>,h){x) = Q[i^{x)Mx),Dh{x),D^h{x),(f>{x),D(f>{x), 
We decompose 

ni m 



ra I 771 \ 

1 - E ^^1) [2 + ^'w] + ^(v- + E ^'20^) ■ (3.11) 
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Then we find 

j=l [ 1=1 

where iV'(s) = F'(w + s) - F'(w) = 0{\s\). It follows that d^Q ^ 0{e^-^) in the 
considered range for the parameters. Now, 

m f m ^ 

3=1 I »=1 J 

since 9hj,w = w' i^z — fj{ey)). Thus dh^Qs, = 0(e'*^'"). We also have dh^Qi = 
0(e'^-^) and 

m 

9D'h,Q = ^D'h.Qi = (l -Y,Q,2)dD^H,S{y) = 0(£^-"), ; = 1,2. 

j=i 

Finally, 

Q = 9*, Qs = Aa^,02 = 0(1), Od^, Q - 2V02 = o(i). 
As a conclusion, using the mean value formula and the facts 

101 _ 02| ^ l^^l _ ^^2| < g-Ht||j0l _ 02|j^^^_^^ 



we readily find the validity of (3.10). In particular, we obtain that for llV'fUoo < 
/3e2-^, ; = 1,2, and 

||Q(V'i,0,h) - Q(7A2,0,h)||p.o < Ce^-ll^i - ^2||^. (3.12) 

Thus, from the contraction mapping principle, we find that for certain /? > large 
and fixed, problem (3.8) has a unique solution = ^(0, h) such that 

||*(<^,h)IU < C{e'' + e^-'Uh.p,^). (3.13) 

The Lipschitz dependence of ^' (3.6) in its arguments follows immediately from 
(3.12) and the fixed point characterization (3.8). □ 

Now, assuming that ||(/'||2,p,a is in the considered range, we substitute ip — "^{(j), h) 
in (3.2) and then obtain that 

rn 

(fl = 1'((/),h) + Q24>j, ij{y, z) = (l)j{y, z - fj{ey)), 

solves problem (3.1) if the vector cj) ~ (01, . . . , 4>j) satisfies the system of equations 
^mJj + F'{v)4>, + (F'(w) + 2)V- + OiA^(«'('/', h) + 0,) + QiS{v) = 0, (3.14) 

in the support of Q2 ■ We want to extend these equations to the entire K.^ x R. We 
recall that in (y, z) coordinates we can write 

Am. - C + - e^zKiey)d, + B, 

where _B is a small operator given by (2.30). It is convenient to rewrite equations 
(3.14) in terms of the functions defined as 

(j)j{y,t) = (j)j{y,t + fj{ey)). 
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We find in coordinates {y,t), 



AmJ, - dl4>j + ^kA, + B]<P, + (3.15) 



where 



- e^K{ey) {t + f,{ey))dtc^ ~ 2e{VKhj{ev) ,VkA4>), 
and, expressed in local coordinates (y,i), y = Yp{j), the operator B'^cj) becomes 



+ e{t + f, (ey)) a], da<t> - 2edih, (ey) d,t4> 



(3.16) 



£^dahj{£Y)dt4> + e^dihjdihji^ej) dtt(t> 



with coefBcicnts evaluated at {ey,et + efj{ey)). The difference between the oper- 
ators Bj and Bj is that the expression for Bj actually makes sense for all {y,t), 
while Bj does only up to |t| < S/e. We set Xo{t) = C(l*l — 10 log e), where, we 
recall, ({t) = 1 for r < 1 and = for r > 2. Then we extend the operator Bj + Bj 
to entire space (y, z) setting 



B, Bj + xoB 



Let us relabel 



and denote 



xUy,t) ~ Qsit + fjiey)) = a\t + hjiey)\-de-s), 



^ji<f>,h){y,t) -.^ 5'(0,h)(y,t + /,(ey)), 



5,(h)(y,t) := Xj35(w)(2;,i + /,), 
(observe that this is the same Sj introduced in (2.53)) 



(3.17) 



:= v{y, t + fj) ^ w{t) + Oj, 



where 9j{y,t) — 0{e ^+^). We have 



dl^i, + A^J, + F'{w,)4,, + 03 B4>j + 03 (i^'(w) - F'iwj))^, 



+ 0i 



(i^'(w) + 2)*(0, h) + N{'f{(f>,h) + 0j) + S'(w) 



in /Cj- X 



where Wj{y,z) = w{z — fj{ey)). Finally, we recast equations (3.18) as 
+ Ak, + F'{w{t))4>j + Sj (h) + Nj- (0, h) = in /Ce X R, 



(3.18) 
(3.19) 
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for all J = 1, . . . , m, where 

%(0,h) = B,(0,) + x,i [ (F'i^j) + 2)vl/,(0,h) + N{^,{^,h) + ] , (3.20) 
with 

B,(0,) = X,3 [Bj(b, + (F'(wj) - F'iw))cj,j^ , given by (3.16). (3.21) 

Wc will concentrate in what follows in solving system (3.19). Wc will do this in two 
steps: 1. solving a projected version of the problem, carrying h as a parameter, and 
2. finding h such that the solution of this projected problem is an actual solution 
of (3.19). We consider then the system, for all j = 1, . . . , m 

dft(bj + ^kA3 = --SjCh) - h) + Cj{y)w'{t) in K., x M, 

(3 22) 

(Pj[y,t)w [t)dt = on ALe , Cj(y) = 



To solve it wc need a suitable invcrtibility theory for the linear operator involved 
in the above equation. Wc do this next. 

4. The auxiliary linear projected problem 
Crucial for later purposes is a solvability theory for the following linear problem: 

40 + A,cA = g{y,t) + c{y)w\t) in JC, x R, 

(4.1) 

<t>[y,t)w'{t)dt^Q for all y G /C^ , c(y) = - ''^ f^' j"" . 

dt 

We have the following result. 

Proposition 4.1. Given p > m and < cr < \/2, there exists a constant C > 
such that for all sufficiently small £ > the following holds. Given g with \\g\\p^a < 
+00, then Problem {^.1) has a unique solution (p with ||0||oo,cr < +oo, which in 
addition satisfies 

\\D^cp\[p,. + P^lloo,. + II0IIOO,. < C\\g\\p^, . (4.2) 

The main fact needed is that the one- variable solution w of (1.4) is nondegenerate 
in L°°(R™) in the sense that the linearized operator 

L{(t>) = /\y(t> + dl(t> + F'{w{t))4), {y, t) e m.^-^ X M, 

is such that the following property holds. 

Lemma 4.1. Let (f> be a bounded, smooth solution of the problem 

L{(j)) = in M.^-'^ X M. (4.3) 
Then <f>{y,t) = Cw'{t) for some C eR. 
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Proof. This fact is by now standard, so we only sketch the proof. The one- 
dimensional operator io(V') = V-"" + F'{w)^p is such that Lo(w') = and w' > 0, 
hence is its least eigenvalue. Using this, it is easy to show that there is a constant 
7 > such that whenever Jj^ ^w' = with £ i/^(M) we have that 

( - F'{w)^' ) dt > 7 / ( l^f + ) dt. (4.4) 



Let (/) be a bounded solution of equation (4.3). Since F'{'w{t)) ~ —2 for all large \t\ 
an application of the maximum principle shows that if < cr < \/2 and io > is 
large then 



|0(y,i)| < C||^||ooe-'^l*l it\t\>to. 
On the other hand, the function 

also satisfies L{4>) ~ and, in addition 



^iy,t) = 4.{y,t) - / ^„'(C)^(j/,C)dC, 



w'{t)(l){y,t)dt = foraU y e M^"\ (4.5) 



Now, the function 



is well defined and smooth. We compute 

Ayifiiy) = 2 Ay4>-4)dt + 2 \Vy4>\'^dt, 



and hence 



^\Ay^~ l_\Vy4>\^dt- l_{\^tf-F'iw)^')dt. 



|2 



(4.6) 



From (4.5) and (4.4), we then get ■^Ayip — -fip > 0. Since ip is bounded, it must 
be zero. In particular this implies that the bounded function 

g{y)= [ wc{OHyX)d(: 

is harmonic and bounded, hence a constant. We conclude that (j){y, t) = Cw'{t), as 
desired. □ 

Proof of Proposition 4.1: We begin by proving a priori estimates. Let < 
cr < \/2. We first claim that there exists a constant C > such that for all small e 
and every solution to Problem (4.1) with ||(/i||oo,iy,o- < +oo and right hand side g 
satisfying < +oo we have 

OO.CT 1 1 0| I OO.CJ < C\\9\\p„^- (4.7) 

To establish this fact, it clearly suffices to consider the case c{y) = 0. By local 
elliptic estimates, it is enough to show that 

II0IIOO,. < q|5llp,.- (4.8) 
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Let US assume by contradiction that (4.8) docs not hold. Then we have sequences 
£ = £n ^ 0, 5„ with ||g„||p,CT 0, 4>n with ||0„||oo,(t = 1 such that 



(4.9) 



dtt(l)n + ^kAu + F' {w{t))(j>,, = g„ in x M, 
(/)„(y, <) w'(i) di = for all y e /C^ . 
Then we can find points (p„, <„) G /C^ x M such that 

We will consider different possibilities. Let us consider the local coordinates for 
/Cg^ around p„ , 

^p„,£„(y) = «^rT^^e„p„(£ny), IyI < ^, 

where l^(y) is given by (2.1). Let us assume first that |i„| < C. Then, the Laplace- 
Beltrami operator of /Ce„ takes locally the form 

^Ke^ = a°-(e„y)9y + e„6"(£„y)5j 

where 

a%{enY) = 5^J + o(l), 6- (£„y) = 0(1). 

Thus 

a%dij<j)n + enb°jdj(j)n + i9tt0„ + F' {w{t))4)n ^ gn{Y,t), |y| < ^, 

where gn{y, t) :— gn (l^(ey), . We observe that this expression is valid for y inside 
the domain £~^lAk which is expanding to entire R^^^. Since (j)n is bounded, and 
(jn in Lf^^{W^), we obtain local uniform W^^'^'-bounds. Hence we may assume, 
passing to a subsequence, that (j)n converges uniformly in the compact subsets of 
to a function ^(y, t) that satisfies 

Ak«-i0 + dtt^ + F'{w{t))4, = . 

Thus (j) is non-zero and bounded. Hence Lemma 4.1 implies that, necessarily, 
(j){y,t) = Cw'{t). On the other hand, we have 

(t>niy,t) w' (t) dt — !► / (i){y,t) w' (t) dt as n — > oo. 
Jm 

Hence, necessarily = 0. But |^„(0,i„)| > ^, and i„ was bounded, the local 
uniform convergence implies 7^ 0. We have reached a contradiction. 

Now, if tn is unbounded, say, t„ — > +00, the situation is similar. The variation 
is that we define now 

0„(y,t) = e"(*"+*V„(y,t„ +t), g^{y,t) = e"(*"+*).g„(y, t„ + t). 

Then 0„ is uniformly bounded, and g„ — > in L^^^{M.^). Now ^„ satisfies 

a%{eny) dij(j)n + dtt(t'n + £nbjiy„+ £ny)dj(pn 

-2adt4>n + F'{w{t + tn) + cr^) 0„ = g„. 
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We fall into the limiting situation 

As,N-i^ + dl4> - 2adt4> - (2 - cr^) = in (4.10) 

with 7^ bounded. The maximum principle implies that = 0. We obtain a 
contradiction that proves the validity of (4.7). 

It remains to prove existence of a solution </> of problem (4.1) with |j(/)||oo,CT < +oo. 
We assume first that g has compact support. For such a g, Problem (4.1) has a 
variational formulation. Let 

■H = i e i?o (/Cj X M) / / (f){y, t) w'{t) dt = for all y e /C^ i . 



H is a closed subspace of Hq {JC^ x R) , hence a Hilbert space when endowed with 
its natural norm, 



+ |Vk.0|' - F'{w{t)) cf" ) AVk, dt . 
is then a weak solution of Problem (4.1) ii (p E H and satisfies 

a(</),^):= / (VkJ-^ic,^ - F'iw{t))(l)^)dVic,dt 
gipdVicdt for aU ipeU. 



It is standard to check that a weak solution of Problem (4.1) is also classical provided 
that g is regular enough. Let us observe that because of the orthogonality condition 
defining Ji we have that 



7/ V d^K, di < a{4>,^) for ah ^eV.. 

Hence the bilinear form a is coercive in Ti, and existence of a unique weak solution 
follows from Riesz's theorem. If g is regular and compactly supported, tp is also 
regular. Local elliptic regularity implies in particular that is bounded. Since for 
some to > 0, the equation satisfied by (j) is 

Acl) + F'{w{t))<j)^ciy)w'{t), \t\>to, y e IC„ (4.11) 

and c{y) is bounded, then enlarging to if necessary, we see that for cr < ^/2, a 
suitable barrier argument shows that |0| < Ce"'^'*', hence |10||p,ct < +oo. From 
(4.7) we obtain that 

<C\\g\\p^.. (4.12) 

Now, for an arbitrary || (7||p,cr < +oo we consider a sequence of compactly supported 
approximations uniformly controlled in || jjp^^ (thus inheriting corresponding control 
on the approximate solutions). Passing to a subsequence if necessary, we obtain 
local convergence to a solution to the full problem which respects the estimate 
(4.2). This concludes the proof. □ 
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5. Solving the nonlinear projected problem 

To solve Problem (3.22) and for the subsequent step of adjusting h so that the 
quantities ci{y) are all identically zero, it is important to keep track of the Lipschitz 
character of the operators involved in this equation. We have the following result. 

Lemma 5.1. There is a constant C > such that for all h' satisfying (3.9) and 
all with 1 1 0' 1 1 2. p. (T < £^"^"7 I = 1,2 we have 

||iV,(</>\hi)-iVj(0^h2)||p,, < C£4-j|hi-h2||w^,p(K)+e2-"||<^i-02||2,p,., (5.1) 

|15,(hi)-5,(h2)|lp,, < Ce3-||hi-h2||w^,p(^). (5.2) 

Proof. We have to check the Lipschitz character of the operators Nj(^, h) in (3.20) 
in the norm || \\p^a. Let us consider each of the terms in formula (3.20). Let us 
consider first the operator Bj(j)j in (3.21). On (j> and h we assume 

\m2,p,a<e^-\ ||h||v^.,p(;c) <M- (5.3) 
This operator has the form in local coordinates 

^3<i>j{y,i) = B{h.j,dikh.j,dih.j,(t)j,Y,t). 

Let us consider the operator Bjipj in (3.16). We sec that the explicit dependence on 
hj comes only from the coefficients anc and &i, more precisely on smooth functions 
of the form a(ey, et + efj{ey)), fj = £,j + hj, so that dhja ~ 0{e). We also find 

dh,X,3 - 0(1), df^.F'i^,) = ^ - ifk - /,)) - Oie'-n. 
Taking these facts into account we then find that for arbitrarily small r > 0, 

and hence 

\\dh,BjC^j\\o,p,. < Ce^-^Uh.p^a . 
Observe that we have as well that 

\\Bj(f>j\\o.p^a < Ce^~^||(/)||2,p,^ . 

Let us consider the dependence on the derivatives of h. We easily check that 

dD^hBj(l)j = 0{e)D^(j,, dD2^hMj = 0{e)D(j)j. 
As a conclusion we find that, emphasizing the dependence on h of the operator Bj, 

||B,(0l,hl) - 8,(02, h2)||o,p < e'-^U' - (t>^\hp.a+S^-^ l|hl - h2||H.2,p(K). (5.4) 

Let us consider the remaining operator in Nj, 

AA(0,h) :-x,i [(F'(w,) + 2)vl/,(0,h)+iV(vl/,(</>,h) + </>,)]. 
We write it as 

AA(0,h)(y,i) =AA(0,V,h,y,i), V = *,(0,h), 
and recall from Lemma 3.1 that II^Alloo = 0(e"'^^). Observe first that 
d^N = Xji [(F'(w,) + 2)+iV'(0 + 0,)] = 0(ee-^l*l), 



d^^Af = xn N'i^ + 0,) - Om + 10,1) = Oie'-^. 
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In addition, we also check that 

(^7V = 0(|V'|e-'^l*l + |0|2) = 0(£4-^). 
Using these estimates, and writing i// = ^E'j h') we find 

iiAA(^\hi)-AA(02,h2)iip^, = ||iv(0^v^h^•)-^(0^^^h^•)llp.^ 

Recalhng now, (3.6), and combining this with estimate (5.4) we arrive to the desired 
result. The proof of (5.1) is concluded. The proof of estimate (5.2) is similar, taking 
into account the explicit form of the error. □ 

Proposition 5.1. Given h satisfying (2.34), problem (3.22) has a unique solution 
4> — $(h) with ||0||2.p.(T < e^""^. Moreover, we have the validity of the Lipschitz 
conditions 

mh') - <i>{h^)h,p,. < C£2-||hi-h2|l,^.„(^). (5.5) 
In addition, we have that 

||Nj(cI>(lii),hi)-N,(a>(h2),h2)||2,p,, <Ce4-n|h'-h'!lH'^-(K)- (5.6) 

Proof. Let T{g) be the operator defined as the solution of (4.1) predicted by Propo- 
sition 4.1. Then we find a solution to (3.22) if we solve the fixed point problem for 

= (01, . . . , <j)N) 

0, =*Bj(0,h) :=-T(5'j(h)+Nj(0,h)) foraU j-l,...,iV. (5.7) 

We will check that the operator 58(0, h) — (*8i(0,h), . . . , Q3jv(0jh)) is a contraction 
mapping in in a ball for the norm || ||2,p,(t- We will do more, checking as well the 
Lipschitz dependence in h. Using the above lemma we conclude that the operator 
*8 is a contraction mapping on the region ||(?!)||2,p,cr < £^~^. Now, using (3.13), 

|S(0)| < x,t((^'(w,) + 2)|*(0)| + |*(0)P + |^,(w)|). 

Thus 

|*B(0)| < e^-^e-'^l*! +Xjie^ + C£2|i:)2.;^(ey)|e--I*l +C£2e-"I*I, 
and hence 

As a conclusion, we can apply the contraction mapping principle, and find a unique 
solution <j> of problem (5.7) such that 

for a suitably large choice of /3. 
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6. The Jacobi-Toda system 

Once problem (3.22) has been solved by = $(h), according to Proposition 5.1, 
the remaining task is to find an h such that for all ^ = 1, . . . , m, we have 

^Ky) = y (-Sflh) +Nj($(h),h)) lu'dt = for all y e JC^. (6.1) 

Using the definition of S( in (3.17), expansion (2.53), Lemma 2.5 and the definition 
of Nj($(h),h), we get 

£-2/,(£-iy) ^h(^A,che + Kiy)fe,) - b2Pe [e-^('''-'''-i) - e-^C^^+i-''^^ " 

+ deih), (6.2) 
where 9i is a small operator: 

for any r > , uniformly on h. The constants bi, 62 are given by 

61= /" w'itfdt, 62= / 6il-w^it))e^'w'it)dt= f 6{l-w\t))e-^'w'{t)dt. 
Jr Jr Js. 



Part I: Recall the relation in (2.37) 



My) = {i 



m + 1 



Pe + he{y). 



Since we want that the functions hi make the quantities Ii as small as possible, it 
is reasonable to find first an h such that the equations, for £ = 1, ... to, 

hi (^A^he + K{y)fe,'j - b2Pe 

be approximately satisfied. We set 



0, (6.3) 



Rt{\i) ■.^<T{AK.hi + K{y)ft 



-\/2{hi-hi_i) _ -\/2(ht^i-hi) 



(6.4) 



where 
and 



a ■.= a,= p^ ^6162 ^ - (log-) ^ , 



R(h) 



i?i(h) 



(6.5) 



.i?,„(h) 

We would like to find a solution h to the system R(h) = 0. To this end, we find 
first a convenient representation of the operator R(h). Let us consider the auxiliary 
variables 

Vl 



_Vrn-l, 

defined in terms of h as 

m 

ve = hi+i - hi with £ = 1, . . . , m - 1, Vm = '^hi, 

1=1 
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with the conventions vq = Vm+i = +00 and define the operators 

- 5i(v) ' 

, S(v) 



S(v 

where we have setted 



S(v) 



_S'm_-i(v) 
5,(v) :=i?,+i(h)-i?,(h) = 



cr (^AK;W£ + A"(y)(pe + W£)) + -| e-^^''+i - 26"^^* + e^^^'-i if 1 < £ 
and 



if ^ = 1, 

< m — 1, 
if £ = m - 1, 



5'm(wm) := ^ i?c(h) = o- f Ak;V,„ + K{y) 



Then the operators R and S are in correspondence through the formula 

S(v) ^ BR(B"iv), 
where B is the constant, invertible N x N matrix 



(6.6) 



B 



-1 1 ••• 
-1 1 ••• 



-1 1 

1 1 1 



(6.7) 



and then the system R(h) = is equivalent to S(v) ~ 0, which setting /3 = 626;^ ^ 
decouples into 



S(v) = a AK:v + A'(y)v + /3K{y) 



+ So(v) = 0, 



(6.8) 



where 



So(v) := -C 



0-V2V1 



-V2t)„ 



c = 



2-10 
-1 2 -1 



4 2 -1 
-1 2 



(6.9) 



(6.10) 



In system (6.8)-(6.9), the second relation and our non-degeneracy assumption force 
v„i = 0. Thus we look for a solution v = (v, 0) of the system, where v satisfies 
(6.8). Rather than finding an exact solution v of S(v) ~ we will find a good 
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approximation. More precisely, by means of a simple iterative procedure, we will 
find for each A; > 1 a function v*^ with the property that 



(6.11) 



Let us find a function with the desired property (6.11) for fc = 1. We consider 
the vector v^(y) defined by the relations 



So(v^) 
Wc compute explicitly 



-PK(y) 



1 



/3 



K(y) {m~£)£ 



, = 1 , . . . , m — 1 , 



3.12) 



and get from (6.8 



This approximation can be improved to any order in powers of cr, as the following 
lemma states. 

Lemma 6.1. Given fc > 1, there exists a function of the form 

where v^{y) is defined by (6.12), = 0, and S,k is smooth on IC x [0,oo), such that 

S(v'^) = 0(a'=) 
as (7 — > 0, uniformly on IC. In particular, 

-k 



:= B 



with B is given by (6.7), solves approximately system (6.3) in the sense that 

R(h'') =0((7'^). 

Proof. In order to find a subsequent improvement of approximation beyond v^, we 
set = + uJi. Let us expand 

S(v^+a;) = a\A,c^^+K{y)uj] + a {A,c^^ +K{yy) + DSoiv^)uj + N{lo), (6.13) 

where 



£)So(v^) V2C 
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2ai —02 
-02 2a3 






■ 




-O.m-3 2a„i_2 ^Im-l 
— 1m-2 2am-l 



(6.14) 



with 



and 



N(v) = |c 



(m-^)£, £ = l,...,m-l, 



(6.15) 



.am(e-^"- -l + \/2< 



The matrix Z?So(v^) is clearly invertible. Let us consider the unique solution 
wi = 0(a) of the linear system 

DSo(v1)l^i = -a{A^v' + Kiy^) = 0{a), (6.16) 
and define — + 1^1. Then from (6.13) we have 

S(v2) = a{AKL^i + K{y)iOi) + N(cji) = ©(a^). (6.17) 
Next we define = + uj2 where CJ2 = 0((7^) is the unique solution of 

- DSo{v^)lo2 = ct(Ak;wi + A'(y)a;i) + N(wi). (6.18) 
Then from (6.13) we get 

S(v3) = cr(AK;W2 + i^(y)<I'2) + N(a;i + 002) - N{uji) ^ 0{<t^). (6.19) 

In general, we define inductively, for fc > 3, v'^^^ = v*^ + Uk where Wfe is the 
unique solution of the linear system 

-Z)So(v^)a;fc = a[Aicujk~i + K{y)ujk-i) + N(a;i H \- ujk-i) 

- N(wi + ---+tjfe_2). (6.20) 
Then clearly ujk = 0(cr'''). Let us estimate the size of S(v*'+-^). From (6.13) we have 



S(v'=+i) = a (Ajcv^ +K{y)v^) + [aiAjc + K) + DSo{v^) + N \ (J^' 
Now, using (6.16), (6.18) and (6.20) we get 
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a{A^ + K) + DSo{v')\ (Y^u:') 

i=l 

k 

+ [(^{^K^i-i + Klo.^i) + i:»So(vi)a;,] 

4=2 

fc 

= -N(wi) - ^[N(wi + • ■ ■ + - N(wi + • • • + a;,_2)] 

1=3 

= -N(wi + ---+a;fc). 



Hence, 



S(v'=+^) cr(AK;a;fe + Kuk) + N(tJi + ■ ■ ■ + uJu-i + i^k) ~ N(a;i + • • • + cjfc-i) 



Finally, the functions •Ci = a-nd 

^^■= a^'^iuj-y^ hwfc.i), fc > 2, 

clearly satisfy the conclusions of the lemma, and the proof is concluded. 



(6.21) 



□ 



Part II: The question now, is how to use the approximation h'" just constructed 
to find an exact h solution to system (6.1). This system takes the form 

R(h) = .9, (6.22) 

where 5 is a small function, actually a small nonlinear operator in h. For the 
moment we will think of g as a fixed function. Since the operator R decouples as in 
(6.6) when expressed in terms of S, it is more convenient to consider the equivalent 
problem 



S(v) = g, 

which, according to expressions (6.8) and (6.9), decouples as 

'1 

S(v) = a AKy + K{y)v + fiK{y) : + So(v) 

1 



3.23) 



(6.24) 



Sm(wm) = o {^KVyn + K{y)v^) = g„ . (6.25) 

Equation (6.25) has a unique solution for any given function gra-, thanks to 
the nondcgencracy assumption. Therefore we will concentrate in solving Problem 
(6.24), for a small given g. Let us write 



where v'^ is the approximation in Lemma 6.1. We express (6.24) in the form 

L<,(w) := -(t[A/cw + A'(?;)J - i:>So(v'^')w = S(v'=) + Ni(u;) + 5, (6.26) 
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where 



Ni(c^) := So(v'^' +Lo)~ So(v^) - -DSo(v'=)w , 



(6.27) 



and So is the operator in (6.10). 

The desired solvabihty theory will be a consequence of a suitable invertibility 
statement for the linear operator L^. Thus we consider the equation 



La-{uj) = g in /C. 



(6.28) 



This operator is vector valued. It is convenient to express it in self-adjoint form 
by replacing the matrix DSo{v'^) with a symmetric one. We recall that we have 




















where the matrix C is given in (6.10). C is symmetric and positive definite. Indeed, 
a straightforward computation yields that its eigenvalues are explicitly given by 



m — 1 



We consider the symmetric, positive definite square root matrix of C and denote 
it by C 2 . Then setting 

w:=C^V, g:=C~ig, 
we see that equation (6.28) becomes 

Laiijj) -.^ -aAjcip ~ A{y,a)ip = g in /C, 
where A is the symmetric matrix 

'e~^'"i 



(6.29) 



A{y,a) = (TK{y)I„^-l + V2C^ 




















C5. (6.30) 



Since 

we have that A is smooth in its variables and 



A(y,0) = -^K{y)C^ 



ai 
as 





O-m-l 



C5. 



(6.31) 



where a£ = £{m — i). In particular, A(y, cr) has uniformly positive eigenvalues 
whenever cr is sufficiently small. 



Our main result concerning uniform solvability of Problem (6.29) is the following. 
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Proposition 6.1. There exists a sequence of values cr = ct^ — > such that 
is invertible. More precisely, for any g G L^{^) there exists a unique solution 
-0 = L^^g G H^(IC) to equation (6.29). This solution satisfies 

<y\\Dli'\\L2(K,) + \mmK) < Ca-^||.g|!i2(^). (6.32) 

Moreover, if p > N — \, there exist C, > such that the solution satisfies besides 
the estimate 

\\DmLP(K)+ \\DK.i'\\L^(K) + \\M\L^(K) < Ccf- \\g\\ lp(^k:) ■ 

In addition, for N ~ 2, we have the existence of positive numbers Vi,V2, ■ . ■ , Vm-i 
such that for all small a with 

\via ~ > C(T~ 2 for all > 1 , i = 1 , . . . , m — 1 , 

for some c > 0, then ip = L^^9 exists and estimate (6.32) holds. 



We postpone the proof of this result in the last section. Assuming its validity, we 
will use it to derive a solvability statement for the Problem (6.22), and to conclude 
the corresponding solvability of system (6.1), hence that of Theorem 1. 

7. Solving system (6.1): Conclusion of the proof of Theorem 1 

7.1. Solving Problem (6.22). Here we refer to objects and notation introduced 
in the previous section. 

Because of the definition of the statement of Proposition 6.1 holds as well 
for the operator in equation (6.28). Choosing a as in the proposition, we write 
this equation as the fixed point problem 

Lo = T{u:) := L'^ (S(v'=) + g + Ni(c.)) . (7.1) 

By construction, we have that 

liS(v'^-)|Up(yc) < Ca\ 
On the other hand, if ||tj||ioo(;i;-) < 5, with 5 sufficiently small we also have that 

l|Ni(a;)|Uoo(^) < C5\ 

and in this region 

||Ni(a;i) - Ni(w2)||l=c(^) < C5\\uji - W2||l~(/c)- 
We observe then that, for v as in Proposition 6.1, 

N — 1 

||r(wi) - T(a;i)||vK2,P(/c) < CSa 5 '^||a;i - w2||lo=(/c), 

and 

\\T{u)U... < Ca-'^-'^ia'^ + S^ + Wgh.^^)). 
Thus if we choose 



fc>2(^^ + j.) 

and g with 

\\9\\lp(k)<ct\ 



V = V{g) 
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then the choice 5 = /icr 2 with ^ large and fixed yields, thanks to the con- 
traction mapping principle, the existence of a unique solution lj to Problem (7.1), 
with 

Let us call uj ~: i^{g). Then, in addition D, satisfies the Lipschitz condition 

N — l 

\\n{gi) -n{g2)\\w^.Pi^n} < Ca^~ " hi - 92\\lp{K) 
for all 51,52 with ^ It follows that the equation (6.23) 

S(v) = g (7.2) 
can be solved under these conditions. In the form 

v*^ + n{g) 

a-\A,c+K)-^9n 

and therefore the equation R(h) — g can be solved for any given small g e L^(/C) 
by means of the correspondence 

S(v) = BR(B-iv). 

This yields the following result. 

Lemma 7.1. Given k > 2{^^-^ + J^); then for all sufficiently small a satisfying the 
statement of Proposition 6.1, and all functions g with 

\\9\\Lp{n) < o-'', 
there exists a solution of the equation 

R(h) = ~g, (7.3) 

of the form 

h = h^- + i7(g), 

where the operator H satisfies 

\\HCgi)\\w-^nn) < 

and 

N — 1 

\\H{gi) - H{g2)\\w^.P(n) < Ca ^ " \\gi - g2\\Lp{K)- 

□ 

7.2. Proof of Theorem 1. We need to prove the existence of h satisfying System 
(6.1). According to expansion (6.2), we have that 

^e-2/,(e-ijy) = R(h) - G(h), (7.4) 
02 

where 

G(h)=-af (?,(y), 
02 

and 9( is the remainder in (6.2). We will estimate this operator. We have that 

0e{y) = i-iy-^e-^ f Ceedh)ie-^y,t)w' (t) dt + e'^ f N,($(h),h)(e-iy,t)«;'dt 
Jm Js. 

" V ' ^ V ' 

Ql(h) Q2(h) 
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where &£ is described in (2.53). Wo have, using Lemma 2.5, 



< Ce-2||C,e,(h)(.,i)||p,„ < Ce'-\ 
And similarly, using Lemma 2.6 we get 

whenever the vector-valued functions hi, ii2 satisfy constraints (2.34). A similar 
argument, using the estimates for the operator Nj{^{h.), h) in Proposition 5.1 yields 

IIQ2(h)|Up(K) < Ce^-\ \\Q2ih') - Q2(h')||Lp(K) < Ce^-llhi - h^\\w^..^,c)- 
As a consequence, the operator G(h) satisfies 

\\G{h)\\L.iK)<Ce'-\ ||G(hi)-G(h2)||i,(^) < Ce^--\\h^-h'\\w-^,^,c), (7.5) 
Thus we need to solve the system 

R(h) = G(h), (7.6) 

which can be rewritten in the form 

R(h^ + q) = G(h^ + q). (7.7) 

We use the operator H{g) defined in Lemma 7.1, and look for a solution of (7.7) 
by solving 

h = H{G{h''+q))=:D{q), (7.8) 
for a sufficiently large fc, in the region 

7^ = {q € W'^P{IC) I \H\w--viK) < /i^'"^"'^}, 

for a sufficiently large fi. From Lemma 7.1 and (7.5), we get that 

||D(qi)-D(q2)|U^.P(;C) < Ga-^-niG(h'= + qi)-G(h'= + q2)|Up(,C) 

' '^^llqi - q2||H'2'!'(/c): 

Hence is a contraction mapping in TZ. Besides we have 

||^(0)|k2,p(o) < Ga'=^^-^ 

I N — 1 

From here it follows the existence of a fixed point q = 0{(t 2 '') for Problem 
(7.8), and hence h = h'^ + q satisfies constraints (2.34) and solves System (6.1). 
This concludes the proof of the theorem. □ 



INTERFACE FOLIATION ON RIEMANNIAN MANIFOLDS 



35 



8. Inverting the linearized Jacobi-Toda operator 

In this section wc will prove Proposition 6.1. The first part of the result holds in 
larger generality. Actually the properties we will use in the matrix function A(y, a) 
are its symmetry, its smooth dependence in its variables on /C x [0, ao), and the fact 
that for certain numbers 7± > 0, we have 

7-|CP<^^A(y,a)^<7+|^|2 for aU ^ e M™-\ (y, a) e /C x [0, ao). (8.1) 



Most of the work in the proof consists in finding the sequence a£ such that 
lies suitably away from the spectrum of L„i, , when this operator is regarded as self- 
adjoint in L^{K). The result will be a consequence of various considerations on the 
asymptotic behavior of the small eigenvalues of as tr — > 0. The general scheme 
below has already been used in related settings, see [19, 20, 21, 16, 17], using the 
theory of smooth and analytic dependence of eigenvalues of families of Fredholm 
operators due to T. Kato[12]. Our proof relies only on elementary considerations 
on the variational characterization of the eigenvalues of and Wcyl's asymptotic 
formula. 

As in the above mentioned works, the assertion holds not only along a sequence, 
but actually for all values of a inside a sequence of disjoint intervals centered at the 

iV-l 

fJiS of width 0{<jf ^ ). The corresponding assertion for N ^ 2 can be made much 
more precise. 

Thus, we consider the eigenvalue problem 

L^(j) = X(j) in /C. (8.2) 

For each cr > the eigenvalues are given by a sequence Aj(f7), characterized by the 
Courant-Fishcr formulas 

Xj{a) ^ sup inf Q^{(j)^(l)) = inf sup Q(t(0, 0), (8.3) 

dim{M)=]-l 06A/^\{O} dirn(M)=j 0gA/\{O} 



where 



We have the validity of the following result. 

Lemma 8.1. There is a number > such that for all < cti < (72 < f * and all 
j > I the following inequalities hold. 

- oi)^ < a^'\,{a2) - CT^'Xjiai) < 2(^2 - ai)^. (8.4) 

In particular, the functions cr G (0, cr*) i— > Aj(cr) are continuous. 

Proof. Let us consider small numbers < cri < cr2 . We observe that for any <p with 
|0p = 1 we have 

o-r^<3'Ti(0,0) - cr2"^QfT2(0,0) = - / '/'^(CTr^A(y,CTi) - (T2"^A(y,cr2))(/) 

JK 



= {(71-^2) [ (j,'^{(J-^A{y,(T)~(T-^d,A{y,a))dp, (8.5) 

JK 
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for some a £ {ai,a2)- From the assumption (8.1) on the matrix A we then find 
that 

From here, and formulas (8.3), incquahty (8.4) follows. □ 

Corollary 8.1. There exists a number 6 > such that for any (72 > and j such 
that 

0-2 + |Aj(cr2)| < S, 
and any cti with i(T2 < Ci < CT2; we have that 

Proof. Let us consider small numbers < cti < (72 such that (7i > Then from 
(8.4) wc find that 

Aj((7l) < Aj((72) + ((72 - (7l) — [Aj((72) - 7 — ], 

(72 (72 

for some 7 > 0. From here the desired result immediately follows. □ 

8.1. Proof of Proposition 6.1, general N. Let us consider the numbers (7£ :— 
2"^ for large i > 1. We will find a sequence of values ae e {(Je+i,(jg) as in the 
statement of the lemma. 
We define 

= {(7 e ((7^+1, (7f) : kerL, ^ {0}}. (8.6) 

If (7 e then for some j we have that Xj{a) = 0. It follows that Aj((7;+i) < 0. 
Indeed, let us assume the opposite. Then, given 6 > 0, the continuity of Xj implies 
the existence of a with < a < a and < Aj((7) < i5. If (5 is chosen as in Corollary 
8.1, and £ is so large that < 5, we obtain a contradiction. 

As a conclusion, we find that for all large £ 

ca.Td{Fi) < N{ae+i), (8.7) 

where N{a) denotes the number of negative eigenvalues of problem (8.2). We 
estimate next this number for small a. Let us consider 0+ > such that 

^^A(2/,(7)e<a+|e|' foraU ^ e M"-\ (y, (7) e /C x [0, (7o), 
and the operator 

L+ = -A^ - — . (8.8) 

(7 

Let A^((7) denote its eigenvalues. From the Courant-Fisher characterization we see 
that A^((7) < Aj((7). Hence N{g) < N^{a), where the latter quantity designates 
the number of negative eigenvalues of L'^ . 

Let us denote by fj,j the eigenvalues of — A^c. Then Weyl's asymptotic formula 
for eigenvalues of the Laplace-Beltrami operator, see for instance [5, 15, 23], asserts 
that for a certain constant Ck: > 0, 

Mi = Cicj'^ + o{jT^) as j +00. (8.9) 
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Using the fact that ^^{o') = l^j ^ ^ ^nd (8.9) we then find that 

N+{a) = C(T-^ +o{a-^) as a ^ 0. (8.10) 
As a conclusion, using (8.7) we find 

card (Fi) < N{ai+i) < Ca^^^^ < C2^— . (8.11) 

Hence there exists an interval (a^, bg) C (ct^+i, a/) such that ae, bi ^ F Ker {L^) = 
{0}, cr e [at, be) and 

0£-a^> jTFT > L.cr^ . (8.12) 

card(r£) 

Let ^ 
We will analyze the spectrum of L^-f ■ If some c > 0, and all j we have 

Af-l 

|A,(af)| >ca,"^, (8.13) 
then we have the validity of the existence assertion and estimate (6.32). Assume 

iV-l 

the opposite, namely that for some j we have |Aj (cr^)! < Sa^ ^ , with 5 arbitrarily 

JV-l 

small. Let us assume first that < Aj(crf) < 5a ^ ^ . Then wc have from Lemma 
8.1, 

Xjiae) < Xj{ae) - {ae-ai) — [Xj{ae)+-f-^]. 

(7e Log 

Hence, (8.12) and (8.13) imply that 

if 5 was chosen a priori sufficiently small. It follows that Aj(cr) must vanish at 
some tr £ (a^, ai), and we have thus reached a contradiction with the choice of the 
interval (a^, bi). 



The case —Sa^ ^ < Xj{ae) < is handled similarly. In that case we get Xj{be) > 
0. The proof of existence and estimate (6.32) is thus complete. 

Let us consider now a number p > — 1 . Now we want to estimate the inverse 
of icTj in Sobolev norms. The equation satisfied by "0 = L~^g has the form 

AKV = 0(a-i)[V' + g] 
for a = at. Then from elliptic estimates we get 

UWy/^.-'iK) < + hhoiic)] (8-14) 

Using this for q ~ 2 and estimate (6.32) we obtain 

WipWw^'^iK.) < Ca^ [\\iP\\l^k) + \\9\\l^k)]<C(t^~~ \\g\\LPiK)- 
From Sobolev's embedding we then find 



H'\\li{ic) < Ca 2 \\g\\LP{ic)- 

I ^ 

N-5 



for any 1 < q < if A^ > 5, and any g > 1 if A^ < 5. If g = p is admissible 



in this range, the estimate follows from (8.14). If not, we apply it for q = ^^z^: 
and then Sobolev's embedding yields 
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\9\\lv{k)- 



for any 1 < s < 



2(jV-l) 



if > 6, and any s > 1 if < 6. Iterating this argument, 
we obtain the desired estimate after a finite number of steps. The proof of the first 
part of the proposition is concluded. □ 

8.2. The case N = 2. Conclusion of the proof. We consider now the problem 
of solving system (6.29) when N — 2. We consider first the problem of solving 



— aA]cip — A{y,0)Tp = g in /C. 

A main observation is the following: the linear system i 
If Ai, . . . , Am-i denote the eigenvalues of the matrix 



(8.15) 

B.15) can be decoupled: 



Q := 



which coincide with those of 



ai 






0-2 







2ai 


-0,2 










2a2 


-as 








-a2 


2a3 






• • • Urn — 3 2fl^j_2 '^ri 

• • • —0,m-2 '2ci„ 

then system (8.15) expressed in coordinates associated to eigenfunctions of Q de- 
couples into m — 1 equations of the form 



- crAicipj 



V2 



/C, j = 1, . . . ,m — 1. 



(8.16) 



When N = 2 this problem reduces to an ODE. K, is then a geodesic of M and 
K(y) will simply be Gauss curvature measured along IC. Using y as arclength 
coordinate, and dropping the index j, Equations (8.16) take the generic form 



-ai;" ~ fiK{y)i} = g in(0,^), 
^(0) = V''(0) - ^'(i), 



.17) 



where fi is given and fixed, and £ is the total length of IC. 

For this problem to be uniquely solvable, we need that fia^^ differs from the 
eigenvalues A = Xj of the problem 



VP(O) = ^i£), ^'(0) = ip'ii). 
More precisely, in such a case we have that the solution of (8.17) satisfies 



IIV^I 



< 



miuj |Aj 



\9\\l^{K)- 



(8.18) 



(8.19) 
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Now, we restate Problem (8.18) using the following Liouville transformation: 



y 



io= VK{y)dy, t = T V^d0, t e [0,^), 

Jo «0 Jo 

^{y)^K{yrK e(i)-^(y)/ q{t) - 



7r^^Hy)Kiy)- 
Equation (8.18) then becomes 

^2 

-e" -q(t)e^ -^Xe in(0,7r), e(0) = e(7r), e'(0) = e'(7r). 
A result in [14] shows that, as j H> oo we have 



•^0 

Hence, if for some c > we have that 



1 47r2j2 



02 
*0 



> ca 2 for all j > 1, 



and a is sufficiently small, then the problem will be solvable, and thanks to (8.19), 
we will have the estimate 

UWl^K) < C<T-5||g||i2(K), (8.20) 

for the unique solution of Problem (8.17). It follows that, under these conditions 
System (8.16) is uniquely solvable and its solution = — (crA^V + ^(y, 0))^^f/ 
satisfies estimate (8.20). 

Now, for a as above, we can write system (6.29) in the fixed point form in L^{IC), 
V' + nV') = -(aAKV^ + A(y,0))-i<?, e L^IC), (8.21) 

where 

T(^) {aA^i; + A(y, 0))-^ [(A(y, <j) ~ A(y, 0)).^ . 

We observe that, as an operator in L^(/C), ||r|| — 0{(j^). Thus, for small a, 
Problem (8.21) is uniquely solvable, and satisfies (8.20). Finally, for the case, 
we argue with the same bootstrap procedure of §8.1. 

The proof of the proposition is complete. □ 
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